Toward a solution of the coincidence problem 

Sergio del Campo*, 1 Ramon HerreraV and Diego Pavon* 2 

1 Instituto de Fisica, Pontificia Universidad Catdlica de Valparaiso, 
Avenida Brasil 2950, Casilla 4059, Valparaiso, Chile 
2 Departamento de Fisica, Facultad de Ciencias, 
Universidad Autonoma de Barcelona, 
08193 Bellaterra (Barcelona), Spain 

Abstract 

The coincidence problem of late cosmic acceleration constitutes a serious riddle with regard to 
our understanding of the evolution of the Universe. Here we argue that this problem may someday 
be solved -or better understood- by expressing the Hubble expansion rate as a function of the ratio 
of densities (dark matter /dark energy) and observationally determining the said rate in terms of 
the redshift. 
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I. INTRODUCTION 



The coincidence problem of late cosmic acceleration, i.e., the fact that the density values 
of dark matter, p m , and dark energy, p^, are of the same order precisely today, constitutes 
a serious challenge to our understanding of the evolution of the Universe. The otherwise 
rather successful vacuum-cold dark matter (ACDM) model offers no explanation. In that 
model, the cosmological constant must be fine-tuned many orders of magnitude for the ratio 
r = p m j p$ to be of order unity nowadays. This partly explains why many authors turned to 
models in which the dark energy density is no longer a constant but evolves with expansion 
-see e.g. [l| for references. 

One way to alleviate the coincidence problem is to provide a mechanism by virtue of 
which r tends to a constant at late times 

aa 

or at least makes it vary slower than the 
scale factor at present time 4J. This can be achieved by introducing an interaction between 
dark matter and dark energy regardless if the latter obeys the dominant energy condition 
or not Q]. 

Most cosmological models assume, for the sake of simplicity, that matter and dark energy 
interact only gravitationally. In the absence of an underlying symmetry that would suppress 
a matter-dark energy coupling (interaction) there is no a priori reason to dismiss it. Further, 
the coupling is not only likely but inevitable [6( and its introduction is not more arbitrary 
than assuming it to vanish. Moreover, it has been forcefully argued that the interaction 
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reveals itself when one applies the Layzer- Irvine equation [7[ to galaxy clusters [8(. Likewise, 
the cross-correlation of galaxy catalogs with cosmic microwave background (CMB) maps 
seems to observationally favor the interaction [9J . Nevertheless, it still remains hypothetical 
and more abundant, varied, and accurate measurements are needed before its existence (or 
non-existence) can be established observationally. 

Cosmological models where dark matter (DM) and dark energy (DE) do not evolve sep- 
arately but interact with each other were first introduced to justify the small value of the 

n 

cosmological constant [10[ and currently there is a growing body of literature on the subject 



-see, e.g. 



ll| and references therein. Further, in some holographic models of dark energy 



the accelerated expansion can be traced to the interaction 12j. Recently, various proposals 
at the fundamental level, including field Lagrangians, have been advanced to account for 
the coupling 
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Obviously, a mechanism that makes r tend to a constant today or decrease its rate 
to a lower value than the scale factor expansion rate ameliorates the coincidence problem 
significantly, but it does not solve it in full. To do so the said mechanism must also achieve 
r ~ We are not aware of any mechanism able to fulfill the latter. 

It is the view of the present authors that, on the one hand, r ought to be understood -at 
least, for the time being- as an input parameter much in the same way as other key obser- 
vational quantities -say, the present value of the cosmic background radiation temperature, 
the cosmological constant H , the age of the Universe to, or the ratio between the number 
of baryons and photons. And, on the other hand, its value must be closely related to Hq 
or, for that matter, to to. Accordingly, if one successfully express r in terms of Hubble rate, 
H, the problem of explaining the r value reduces to the problem of explaining the present 
value of H. This is interesting because while we are still lacking model-independent data of 
r at different redshifts, two sets of observational data values -albeit scarce- of H at various 
redshifts are now available. The first data set was obtained by computing differential ages 
of passively evolving galaxies in the redshift range 0.1 < z < 1.8 14J. The other data set 
is based on measurements of 192 supernovae type la (SN la) and 30 radio galaxies up to 



redshift 1.2 
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Before proceeding, we wish to emphasize that, to the best of our knowledge, no model- 
independent data of the evolution of the energy densities of matter and dark energy do exist. 
To obtain such evolution (e.g., with respect to redshift) one must specify some cosmological 
model. It is to say, the field equations together with the equation of state of dark energy 
and some assumption about whether the components interact or not. If they are assumed 
to interact, one must also provide an expression for the interaction. 



II. BASIC RELATIONS 



Let us consider a spatially flat Friedmann-Lemaitre-Robertson- Walker universe domi- 
nated by a two-component system, namely, pressureless dark matter and dark energy, such 
that the components do not conserve separately but interact with each other in a manner 
to be specified below. The energy density and pressure of the dark energy, assuming it is a 
quintessence field, are given by p^ = \4> 2 + V (</>), and = \ <fi 2 — V (</>), respectively. If the 
dark energy is a phantom field we have instead, = — ^<p 2 + V(4>) , and P^ = — ~ 4> 2 — V(4>). 
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The upper-dot stands for derivative with respect to the cosmic time and V(4>) denotes both 
the quintessence field potential and phantom potential. As is usually done, we postulate that 
the dark energy component (either quintessence or phantom) obeys a barotropic equation 
of state, i.e., = w p^ with w a negative constant of order unity and lower than —1/3 (a 
distinguishing feature of dark energy fields is a high negative pressure). 

We assume that the dark matter and dark energy components are coupled through a 
source (loss) term (say, Q) whence the energy balances take the form 



In what follows we shall consider Q > 0. On the one hand, this choice ensures that r de- 
creases monotonously with cosmic expansion (as expected in usual models of cosmic struc- 
ture formation) and that around present time it varies very slowly, i.e., slower than the scale 
'actor. On the other hand, it shows compatibility with the second law of thermodynamics 
1^ | . We note, parenthetically, that if Q were negative (i.e., if dark matter decayed into dark 
energy), it would exacerbate the coincidence problem since the ratio r would decrease faster 
than in the ACDM model -see Eq. ([3]) below. 

In virtue of the above expressions and Friedmann's equation, 



p m + 3Hp 



m 



and 




(1) 



3H 2 = K 2 (p m + p ) 



(k 2 = 8tiG), 



(2) 



the time evolution of r can be written as 



r = 3Hr w + 



K 2 Q(r + lf 
g H 3 r 



(3) 



Using the relationship f = H dr / dH , where 



H 



(Pm + P<t> + P<t>) 



3 1 + w + r 



H 2 , 



2 1 + r 



Eq. ([3]) becomes 



Here 



9 = -2r- (1 + r) 



[1 + w + r) 



w + 



K 2 Q{r + lf 
9H 3 r 



(5) 



Clearly, Eq. can be integrated whenever an expression for Q in terms of H and r is given. 

Notice that to have dr/dH > (a very reasonable assumption since one expects that 
H(t < to) > Hq) the relationship 

\w\ > 6 

9 H r K ' 

must be fulfilled. 

We note in passing that for the ACDM model (i.e., w = — 1 and Q = 0) one follows 
H = H [(l + r)/(l + r )]V 2 . 

III. TWO CHOICES FOR Q 

In the following, we focus on two phenomenological expressions for Q previously consid- 
ered in the literature. As we have carefully checked, both models fit very well the observa- 



Q-fl, E Q- As depicted in Figfl], it leads to a 



tional data of Daly et al. obtained from selected sets of SN la and 30 radio galaxies [15 ]. 
In particular, the maximum variation of does not exceed 23% in the redshift interval 
< z < 1, thus showing compatibility with the data analysis of Daly et al. [ijj]. 
(i) Q = 3a H (p^ + p m ) with a a positive-definite, dimensionless constant. 

This choice has been studied in detail 
constant but unstable ratio r at early times (high redshifts) and a lower, also constant but 
stable (attractor) ratio at late times. When the dark energy is of phantom type (not shown 
in the figure), the evolution of r is also similar. Further, it fits rather well data from SN 
la, cosmic microwave background radiation (Wilkinson Microwave Anisotropy Probe 3 yr), 
and structure formation -namely, Sloan Digital Sky Survey, and Two Degree Field Redshift 
Survey- provided a < 2.3 x 10~ 3 [171 ] . Likewise, it shows compatibility with cross correlations 
of galaxy catalogs with CMB maps provided a < 0.1 [9j. In this case Eq.([5]) reduces to 
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FIG. 1: Evolution of the ratio r = p m / P<f> with redshift for different values of the parameter a. In 
drawing the curves we have fixed tq = 3/7 and w = —0.9. 
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-It 



1 + r 



w + a 



[r + If 



(7) 



1 + w + r 

where we used Eq. (j2J). The condition given by Eq. (jBj) leads to restriction \w\ > a (l+r) 2 /r 
Using last equation in (jlj) and integrating, the Hubble function is found to be 



H = H exp [/(r)-J(ro)], 



(8) 



where J(r) stands for the real part of J(r) with 



Kr) 



— ln(l+r) — — ln[— w r — a(l +r) 2 ] 



1 



w + 2 



tan 



2a(l + r) +w 



■ (9) 



2 w(4a + w) a/— w (4a w) 

Figure [2] shows the dependence of the Hubble function on the densities ratio for a = 
10~ 4 and two values of the equation of state parameter, w. Likewise, we have plotted the 
prediction of the ACDM model. In all the cases tq = 3/7. For other a values compatible 
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with the restriction a < 2.3 x 10 3 set by the Wilkinson Microwave Anisotropy Probe 3 yr 
experiment the corresponding graphs (not shown) are rather similar. 
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FIG. 2: Evolution H vs the ratio r — p m /p^ as given by Eq.© with a = 1CT 4 for w = -0.9 
(quintessence, top line) and w = —1.1 (phantom, bottom line). Also shown is the prediction of the 
ACDM model (middle line). In drawing all the curves we have fixed ro = 3/7. 

Unfortunately, as said above, no model-independent data of r at different values of H (or 
z) exists whereby, for the time being, we cannot directly contrast this model with observa- 
tion. This is why we turn to determine the dependence of the Hubble factor with redshift 
and compare it with the two available observational data sets H vs z. The first model- 
independent data set was obtained from the study of the differential ages of 32 -carefully 
selected- passively evolving galaxies in the redshift range 0.1 < z < 1.8 [ijj]. The age of 
each galaxy was found by constraining the age of its older stars with the use of synthetic 
stellar population models. The differential ages roughly yields dz/dt, then H(z) is given by 
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H = —dz/[(l + z) dt] -see right panel of Fig. 1 in Ref. [14J]. The second model-independent 
data set (see lower panels in Fig. 9 of Daly et al. [lsj]) was obtained by applying the model- 



independent analysis method of Ref. [18 1 to the coordinate distances of 192 SN la of Davis 



20 1 



et al. [l9j and 30 radio galaxies of Ref. 

To obtain the expression for H(z) of this model we combine Eq. 
([3]) in terms of z which is 



with the integral of 



r(z) 



Re< — 

2a 



-w + 7 tan ( A — — 7 ln(l + z) 



(10) 



where Re specifies the real part of the corresponding quantity, A = tan -1 [(w+2a(l+r ))/7] , 
and 7 = a/— w (w + 4 a). 

Figures [3] and H] show the dependence of the Hubble expansion rate on redshift predicted 
by the model for a = 10~ 4 and two values of the equation of state parameter, w. For 
comparison, we have also plotted the prediction of the ACDM model. The observational 
data in Fig. [3] are taken from Simon, Verde, and Jimenez [jjj] (full circles) and Daly et al. 
(full diamonds), and in Fig. H]from Daly et al. (cf. lower panels of figure 9 in Ref. 15]). 
As can be seen in both figures, the fit to data does not vary significantly between graphs. 
Given the scarcity of data it is not worthwhile to compare the \ 2 °f the different curves. 
(ii) Q = 3 (3 H p m with ft a dimensionless positive-definite small constant. 

Upon this choice, considered in Ref. 4J, the ratio r does not tend to a constant (i.e., 
unlike the previous case no attractor exists 2l|) but, as Fig. [5] shows, it varies very slowly 
at late times -by very slowly we mean that | (r/r) H , "soft coincidence"- whereby the 
coincidence problem gets also significantly alleviated. 



Likewise, in this instance, Eq.(jSJ) reduces to 
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-2r 



1 +r 



1 + w + r 



[w + ft{r+ 1)] 



(11) 



where we used Eq.© to obtain Q = 9/3 r H 3 /[k 2 (1 + r)]. Now, the condition set by Eq.flHJ) 
boils down to \w\ > ft (1 + r). 

With the help of Eqs.(pE]) and ffTTj) the Hubble function can be cast as 
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FIG. 3: Evolution H vs z with a = 10 -4 for w = —0.9 (quintessence) and w = —1.1 (phantom). 
Also shown is the prediction of the ACDM model (solid line). In all the cases we have fixed 
ro = 3/7 and Hq = 71km/s/Mpc. The data points with their la error bars are borrowed from 



Simon, Verde, and Jimenez, Ref. 14| (full circles) and table 2 of Ref. 



151 ] (full diamonds). 



H=H,j£y (12) 



where 
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FIG. 4: Same as Fig. [3l The upper panel shows the combined sample of 132 SN la of Ref. [19(] and 
30 radio galaxies of Ref. 20]. The bottom panel shows the 30 radio galaxies only. The data points 
with their la error intervals and the best-fit curve (big solid line) in both panels are borrowed from 



Daly et al, Ref. 



15] 



J(r) = v / + r^" 1 ^ 2 ^ [/3(1 + r) + w f-MMf>+»)] 



(13) 



Numerically, the dependence of the Hubble expansion rate upon r is very close to the 
previous case (i) for (3 values similar to a. This is why we do not show it here. 

The next step is to compute the dependence of the Hubble expansion rate on redshift. 
This can be done by combining Eq. ()12p with 



r(z) 



£r 



(3r -(l + z) 



-3? 



1 + r ) + w\ 



(14) 



where £ = —w — (3 > 0. The latter expression follows from integrating Eq. in terms of z. 

Again, the corresponding figures for H(z) is rather similar to Figs. [3] and |4] whereby we 
do not depict them here. 
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FIG. 5: Evolution of the ratio r = p m j P& with redshift for different values of the parameter (5. In 
drawing both curves we have fixed ro = 3/7 and w = — 0.9. As it is apparent, r/r ~ at late 
times. 

IV. DISCUSSION 

The current value of the ratio r = p„J is of order unity (more precisely, about 3/7 -see 
e.g. 22]). This represents the "coincidence problem" that the ACDM model seems unable 
to account for. Nowadays, r cannot be derived from first principles and like some other key 
cosmic quantities it has to be considered an input parameter. It seems obvious, nevertheless, 
that in Einstein relativity it must be closely tied to the present value of the Hubble expansion 
rate, Hq, another quantity whose value we cannot explain but that contrary to ro it does not 
mean any puzzle. Therefore, a possible avenue to someday solve the coincidence problem 
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may well be to express r as a function of H (or vice versa) in models such that r at late 
times either tends to constant or varies slowly enough (i.e., slower than in the ACDM model). 
Thus, the problem of explaining r reduces to the less acute problem of explaining H . 

In this paper we have taken a first step in that direction. Specifically, we have considered 
two phenomenological models, previously studied in the literature, in which dark energy and 
dark matter do not evolve separately but interact with each other in such a way that either r 
stays constant at late times (case (i)) or varies very slowly (case (ii)) and, in both instances, 
we have related H to r. Since mo del- independent data relating these two quantities are not 
currently available we have expressed r in terms of z [Eqs. ( ITU]) and ( |T4|) . respectively] and 
resorted to the observational data of Refs. [14| and 15j that link H to z. Regrettably, given 
the scarcity and dispersion of the data (cf. Figs. [3] and H|) they are unable to discriminate 
between models . 

Nevertheless, it is fair to say that in order not to spoil the well established primeval 
nucleosynthesis scenario, fl^ should not exceed 5% at that early epoch -see Bean, Hansen, 
and Melchiorri j^ . At any rate, as Figs. Q] and [5] show, the models considered in this 
work are consistent with this constraint which, in any case, was determined using specific 
non-interacting dark energy models. 

It is to be hoped that in the not far future abundant and accurate observational data, in 
an extended redshift range, may help determine H in terms of z and solve (or, at least, help 
to better understand) the coincidence problem of late acceleration. Yet, it may well be that 
the ACDM model comes to fit the data better than any interacting model. In such a case, 
we should conclude that the coincidence problem is just that: a coincidence and, therefore, 
no longer a problem. 
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